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Abstract 

We extend earlier constructions of the effective action for neutrino-photon scattering, using 
the connection between low-energy neutrino-photon and photon-photon scattering together 
with the known effective Lagrangian describing low-energy photon scattering in QED. We 
use this effective action to calculate analytic expressions for the low-energy cross section 
for the (unpolarised) processes vv ^ 57, ^7 — > z^ + 47 and 77 — > uu + dj. As a byproduct 
we derive compact expressions for the A-body phase-space integrals for massless particles, 
including those having non-trivial tensor-structure. 

1. Introduction 

Neutrinos and photons may be, with gravitons, the only particles which are massless, 
or very nearly so. As such, they are the only degrees of freedom which arise at extremely 
low energies within the vacuum sector (and possibly within other sectors) of the Standard 
Model (SM). This makes the study of their low-energy interactions a theoretical labo- 



ratory for very-low-energy Standard-Model physics. This study may also have practical 
applications, despite the extremely weak strength of the interactions, because neutrinos 
play a unique role within the extreme environments found in astrophysics and cosmology. 
Precisely because of their weak couplings they are often the mediators of dynamically in- 
teresting processes, for instance by being responsible for heat and momentum transfer, 
especially in the late stages of stellar collapse. 

2^2 processes — like vv — > 77 and z/7 — > z/7 — were first studied long ago and 
were found to be highly suppressed [1],[2]. The suppression arises because Yang's theorem 
[3] prohibits the coupling of two photons to a state of angular momentum one, and this 
ensures that the 0(l/M^) contribution to the amplitude for these 2^2 process must 
be zero (but see [8]). The dominant contribution therefore arises at 0(1/Af£), making 
it smaller than the nominally negligible processes which arise at higher order in a, such 
as z/z/ — > 37 and z/7 — > v + 27 [4] . This observation has stimulated more detailed studies 
of these reactions (both at low energies and at energies above the electron mass [5], [6]), 
as well as generating searches for practical applications of these 2 — > 3 processes, such as 
within stars or in neutrino- photon scattering in the presence of magnetic fields [7]. 

As pointed out in ref [4], the effective action for these 2 — > 3 processes is related 
to the known Euler-Heisenberg effective action for four-photon scattering through the 
replacement of one of the electromagnetic tensors, F a p, by a neutrino 'field strength', of 
the form N a p = d a (v^p 7 L z^ — (a <-> /?). In ref. [9], this connection was examined in some 
detail: a particular combination of Feynman diagrams was employed to explicitly show the 
mapping between the calculation of the Euler-Heisenberg action and that governing the 
interactions of neutrino-antineutrino pairs with three photons. 

In this note we have two goals. Our main new result is to extend this treatment 
to the next least complicated case: that of 2 neutrinos interacting with 5 photons. We 
do so by computing the relevant terms of the low-energy effective neutrino-photon la- 
grangian, and use them to calculate analytic expressions for the low-energy neutrino- 
photon scattering cross sections. In the centre-of-mass (CM) we find these to be of order 
a ~ a 3 G% i? 18 /(27r) 6 mg 6 , where E denotes the CM scattering energies, as compared with 
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the 2 — > 2 result: a ~ 200 aG^E 6 /n 3 . In principle, the numerical factors are such that the 
2^5 processes can dominate the 2 — > 2 processes for E ~ m e (which, of course, lies at 
the limit of validity of the low-energy approximation), although we know of no practical 
application of this observation. 

Our secondary goal is not so much new as it is explanatory. Preparatory to describing 
the above results we recast the argument for the suppression of the 2 — > 2 processes into 
a more modern effective-lagrangian language. We also rederive the connection between 
the electromagnetic and neutrino scattering processes within this context. Although these 
are old results, we hope that their recasting in this way may suggest more applications 
elsewhere. 

Our presentation is as follows. In the next section, §2, we review the low-energy 
limit of neutrino-photon scattering, rederiving both the connection to the Euler-Heisenberg 
effective lagrangian and the suppression of 2 — > 2 processes. This is followed in §3 by the 
derivation of the low-energy effective action for 2 — > 5 neutrino-photon interactions. §4 
then applies this action to compute the three 2^5 cross sections: vV — > 57, 1/7 — > v + 47 
and 77 — > uV + 37. We conclude in the last section, §5, with comments and final remarks. 
Finally, an appendix describes an efficient method for computing the relevant iV-body 
phase-space integrals which are encountered. 

2. Low-Energy Neutrino-Photon Scattering Revisited 

Because the SM contains no direct (tree-level) couplings between neutrinos and pho- 
tons, the starting point for calculating their very-low-energy 1 interactions is the SM 
description of their couplings to other particles. These other particles then generate the 
effective neutrino-photon interactions once they are integrated out to produce the very- 
low-energy theory. The effective couplings are nonrenormalizable, in the sense that they 
are proportional to inverse powers of the masses of the particles which were integrated 
out to obtain them. Our main interest in what follows is in the dominant interactions at 

1 We use the name 'very-low-energy' to mean energies below the electron mass, in order to distinguish this 
from other potential notions of 'low energy', such as E<tiM w . 
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very low energy and so we focus on integrating out the lightest particles. Since the two 
lightest particles which couple to both neutrinos and photons are electrons and muons, we 
concentrate our attention on these. 



2.1) The Weak-Scale Effective Theory 

At energies below the VF-boson mass, the couplings of neutrinos and photons to 
charged leptons are described by the effective lagrangian obtained from the SM by in- 
tegrating out the top quark and the electroweak gauge bosons, W and Z. The resulting 
effective interactions which are of most interest in what follows are those which are sup- 
pressed by the fewest powers of M w or M z . Those involving just neutrinos, photons and 
charged leptons, obtained by matching to the SM at ji = M w , are given by: 2 



£ wk (^ = M w ) = e J£ m + -| J2 ( w klvlLVi) Ll lmn + O U-\ , 

V klmn ^ w ' 



(1) 



where k, l,m,n = e, fx, r run over the three lepton flavours, and the charged-lepton currents 
are given by 

^klmn 1^ {^klmn ~\~ ^klmnlb) 

k 

The effective couplings, Vkimn and dumm as found from tree- level matching are given by: 



Vklmn(H = M w ) = SknSlm + Ski S mn ( ~ 7i + 2s w 
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and a k i mn (/i = M w ) = S kn S im - -Ski S r , 



(3) 



at tree level, where s w = sin 9 W is the sine of the weak mixing angle. 

Radiative corrections are easily incorporated into this language. Those loops involving 
high energy degrees of freedom (those involving particles having masses as large as M w or 
larger) are included by matching to the SM with higher-loop accuracy. Intermediate-scale 



2 Like all God-fearing people, our conventions are: ?7 M „ = (— ,+,+,+) and ip = ilft^ 7° . 

4 



loops (involving particles having masses between m e and M w ) are obtained by running 
the effective theory down to each new particle threshold, and then matching across this 
threshold. 

One such high-energy loop generates an effective coupling between neutrinos and 
photons which is proportional to [4]: 



4a (G F 



As we shall see, this particular higher-dimension interaction is not generated when lighter 
particles are integrated out, and so it is the dominant contribution to low-energy 2^2 
photon- neutrino scattering even though it is suppressed by four powers of M w ? 

Imagine now writing down the effective theory at scale ii = m e , just before integrating 
out the electron. The only particles in this low-energy theory are the electron, photon 
and neutrinos. Conservation of electric charge and lepton numbers require the lowest- 
dimension neutrino couplings to electrons in this theory to again have the form of eq. (1), 
although now restricted to electrons and neutrinos. Furthermore, since all of the neutrino 
interactions in this effective theory must vanish in the limit where the W and Z become 
infinitely massive, they must be proportional to at least one factor of G F . 

It follows that the dominant low-energy interactions in the effective theory at this 
scale can differ from the electron terms of eq. (1) only through the values taken by the 
the coefficients Vkiee{n = Tn e ) and a,kiee{n = m e ). Happily enough, it also happens that 
Vkieein = "7-e) cannot differ from its value, eq. (3), at ii = M w , because the current e"f a e 
is conserved, and so does not get renormalized. 

At low energies the sole contribution of physics between \i = M w and \i = m e therefore 
is to the running of the coupling ci^iee 

(and of the electric charge, e) between these scales, 

to all orders in all other SM couplings. 



In very-low-energy scattering applications it is the effective couplings renormalized at the electron mass 
which are required, so /i=m e is used in eq. (4). 
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2.2) Matching at m e 



Next integrate out the electron itself to obtain the effective theory of photons and 
neutrinos only. At this stage we keep effective interactions having more than the minimal 
dimension, because these receive their largest coefficients when the lightest possible particle 
- the electron - is integrated out. We obtain in this way all contributions to low-energy 
neutrino-photon physics which are (3(1/M^mf), for all p. 

To this order the new contributions to the effective neutrino-photon interaction la- 
grangian obtained by matching across the electron mass threshold is therefore given by 



G 

£elth(> = TO e ) = ^ {wklvlLVl) (vklee (e ^ e) + Clklee (e 7^^) ) , (5) 



where (X^) represents the expectation of the operator X M , obtained by integrating out 
the electrons, weighted by the QED lagrangian: 4 



(X») = J VeVeX»(e,e) exp i J d^x (£ kin - ie A» e Y e) 



(6) 



Eqs. (5) and (6) contain the nub of the main results, because it permits the following 
two conclusions: 

• Suppression of 2 — > 2 Processes: 

As is easy to show, all operators involving only V^^ L v - as opposed to z77 Q 7 L 9^ v - 
and two electromagnetic fields vanish on using the equations of motion for the neutrino and 
photon fields, and so are redundant in the sense that they may be removed by performing 
a field redefinition. The only possible lowest dimension operator for 2^2 processes 
(dimension 6) turns out to be of the form of eq. (4). 



4 A notational aside is in order here, since eq. (5) gives the impression that (X M ) does not involve an 
integration over the electromagnetic field as well as the electron field. In reality this expectation denotes 
the usual matching procedure: the difference between the average calculated with electrons and photons in 
the theory just above m e , and the average calculated with photons only in the effective theory just below 
m e . This distinction plays no role in the present discussion. 
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It remains to show that operators of this form are always suppressed by at least two 
powers of G F . We have just argued that the right- hand-side of eq. (5) is explicitly propor- 
tional to the neutrino current, VklnlL^i^ and so cannot contribute to an operator with a 
derivative embedded within the neutrino bilinear. This is why integrating out the electron 
does not generate the operator, eq. (4), with a coefficient proportional to G F /m 2 e . The 
same argument also precludes generating such a term when the other charged leptons are 
integrated out. Charged-current interactions of neutrinos with quarks, on the other hand, 
can be linear in u, and so need not be proportional to VklnlL^i- Nonetheless, conser- 
vation of quark flavour only permits these interactions to contribute to neutrino/photon 
scattering at second order in G F . 

• Connection with Photon-Photon Scattering: 

Since the electromagnetic interactions preserve parity (V) and charge conjugation (C), 
these symmetries may be used to further organize the contributions to £ e ith- In particular, 
these symmetries imply that any term in £ e ith involving an odd power of receives 
contributions only from the vector current, (e7 M e), while those involving even powers of 
F^ u arise purely from the axial current, (£7^75 e). 

Furthermore, since the vector current, e7 M e, is also the electromagnetic current for the 
electron effective theory, its expectation may be expressed in terms of the Euler-Heisenberg 
effective lagrangian, W EH [A] , for photon-photon scattering below m e [4] , [9] : 



e \8A^ 



r r r , ^ (7) 

where Z[A] = e lWEH ^ = / VeVe exp % / d 4 x (C kin - ieA^e^ej 



For instance, since the quartic contribution to the Euler-Heisenberg interaction is 
given by: [10] 

C { E i = -^-j [5(F F F V ) 2 - UF^F» x F Xp F^] , (8) 

it follows that the dominant contribution to £ e ith involving two neutrinos and three elec- 
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tromagnetic fields must be [4], [9]: 



(3) = e(\+2sl)a (G 
elth ' 90nm 4 e V\/2. 



5 (iV F^) (F\ p F Xp ) - 14 {N„ v F vX F Xp F™)] , (9) 



with N a p = d a (vypiL") -(«<->/?). 

This method clearly works in general: to obtain any term involving an odd power of 
F pv in £ e ith, replace one power of e by G F (| + 2s^) /y/2, and sum all possible ways of 
replacing one electromagnetic field strength by N^. 

2.3) The 2-^5 Effective Lagrangian 

The next simplest neutrino-photon interaction which is related in this way to the 
Euler-Heisenberg action describes 2^5 processes, like Vv — > 57. It is related to the sixth 
order term of the EH action, which is given by [11]: 



,(6) 
'EM 



net 



315 m 8 e 



[9(F Q/3 F^) 3 - 26F pu F uX F\ p F pp (F Q/3 F^)] 



(10) 



Now, given our previous arguments we can read off the effective two-neutrino/five- 
photon operators directly. We find: 



eff " 315 y/2 ™e 



6 • %F aP F^) 2 F^N^ 



4 • 26F^F» x F Xp N p » (F^F^) - 2 ■ 26F^F» x F Xp F p » {F af3 N af3 ) 



(11) 



A similar method will also give the even powers of F^ u in £ e ith given the expression 
for the axial- vector /vector current correlation in QED 5 . 



After completing this paper it was brought to our attention that, in ref. [13], the expression for (e7 M 7 5 e) 
has been worked out, up to fourth order in the fields. We thank H. Gies for pointing this reference out to 
us. 
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3. Cross sections for v v — > 5 7 and crossed processes 



We next apply the lagrangian, eq. (11), to compute the 2 — > 5 processes 17z/ — > 57, 
1/7 — > z/ + 47 and 77 — > 77z/ + 37. This is a straightforward, if tedious, exercise within 
the effective theory, requiring only the Born approximation using interaction (11). (This 
should be contrasted with the difficulty of extracting the low-energy limit of the scattering 
amplitude, computed directly from the higher-loop graphs involving the weak interaction, 
eq. (1), and QED. As is usually the case with effective lagrangians, the payoff in simplicity 
is much larger for nonleading contributions.) In performing this calculation we employed 
the symbolic manipulation program FORM [14], which reduced the squared amplitude to 
its final form in under ten minutes on a desktop PC. We briefly sketch the method of 
computation below. 

We require, then, the matrix elements of the effective interaction, eq. (11), which we 
write as follows: 

Cl~^ = gN^T^>->^d ai A Pl ■ --d a5 A P5 

where g is the factor premultiplying the square bracket in eq. (11) and T^ 1, "' ,a ^ 6 repre- 
sents the polynomial of momenta in the effective interaction. In terms of these quantities 
the matrix element relevant to Vv — > 57, for instance, becomes: 



(12) 



where = (0\N^ u \uu) and 

yai/3i," -,as/3s _ ^ y ^ a,r i ''"' a7r 5 (13) 
7r€S 5 (l,",5) 

is the permutation-summed tensor contracting the fields together. The e^'s are, as usual, 
the photon polarisation vectors. 

After squaring and doing the spin sums, the following phase-space integral is required 
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to obtain the total cross section: 



~W ' ' '^t k?lkJ1 " ' ^ kl ~ w) ' (14) 



A general technique for evaluating integrals of this form is given in the Appendix. 
• vv ^ 57; Using the integrals of the Appendix gives the final result for vv — > 57: 

~ 6.87 • 10" 38 I barn, 
\m 2 J 



(15) 



where s = —(p + p) 2 = —2p ■ p is the usual Mandelstam variable, equal to s = AE 2 in the 
centre-of-mass frame. 

• z/7 — > z/ + 47; A similar exercise, after crossing the external lines, yields 

a(v-f — > z/ + 47) = a(z/7 — > 17 4- 47) 

13- 163-2339 5 2 ^ 2 fl n 2 \ 2 / s > 9 



(2tt) 6 2 8 3 10 5 5 7 4 11 e F V2 / W/ ( 16 ) 

:-: 8.68 ■ 10" :!,s (-^1 barn. 



77 — > vv + 37; Crossing the other neutrino yields, 



o x 797549 5 2^2 A n 2 1 / •• 

a(77 -> ^ + 3 7 ) = (2?r)6 27395574n <* m e G F ^- + 2, w j (-^ 

~ 8.38 • 1(T 38 ( -At ] barn. 



(17) 



4. Conclusion 



We have constructed the effective interaction which governs the interactions of five 
photons and two neutrinos using the general connection between the effective action for 
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neutrino-photon interactions at lowest order in G F and the known Euler-Heisenberg effec- 
tive interaction for photon-photon scattering. As an application we have computed the 
two-body cross sections whose low-energy limits are given in terms of this effective in- 
teraction. While these cross sections are likely to be too small to be of any astrophysical 
relevance, it is interesting that the effective interaction for such a high-order process can be 
obtained with such minimal effort. It is also noteworthy that a seventh-order process can 
compete with the two-body scattering Vv — > 27. Finally, our expressions were obtained 
by evaluating multi-body phase space integrals, for which we have presented an efficient 
method of computation. 
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Appendix I. Phase Space Integrals 

Our goal in this appendix is to describe how to evaluate the integral, eq. (14), which 
we reproduce once more for convenience: 



r rfih, rl s k 



We will find, through the integral representation of the delta function, that we are able to 
reduce the problem to one of taking derivatives of a suitable (simple) integral. Since the 
general expressions are extremely lengthy, and particular results are not difficult to obtain 
with the help of a symbol manipulation program once the prescription is known, we will 
only provide a detailed recipe for evaluating these integrals. 

We proceed by using the Fourier representation of the delta function to factorise this 
integral into products of terms having the form 



ik-x 
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so that 
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2«i...a ro ;7i-7m( 1i ,) = y AjL J«i7i ( x ) . . . J*mjm ( x ) e -*«»-*_ ( 18 ) 

The integral defining J a ^ is easily performed as follows: 
where we use the integral 



(19) 



(20) 

We have ensured the convergence of this integral through the appropriate e prescription, 
taking x 2 = —(x ) 2 + x 2 + i sgn(x°)e. The e term in x 2 forces the incoming momentum 
to be future-pointing, and allows the x integrals to be done unambiguously. Notice that 
our result for J al3 is traceless, as is required, since the ks are null. 

With J af3 in hand, expand the integrand of (18) to obtain a sum of integrals of the 

form 

r „—iw-x 



(Txx ai ---x^ 



(x 2 )™' 

where w is a future-pointed, timelike four-vector. These integrals can all be evaluated by 
differentiating 

/g— iw-x 
d 4 x (~2yn , (21) 

with respect to w, so that finding this integral reduces the problem of calculating (14) to 
one of expanding a polynomial and taking derivatives. 

Let us evaluate (21). If we explicitly factor the denominator, go to the rest frame of 
w, and consider the t integral first, we need to consider 



1 1,2(1 1 



oo poo „+iut 



l m (w) = / eTfi / drr 2 / dt 



e 



DC 



'[-(t-ie) + r] m [(t-ie)+r] 



i p poo poo „+iu)t 

- / d 2 n / dr r 2 / dt 



27 ' J-oo ' i-oo (-)"»[ *-(r + ie)]™[*+(r-ie) 

/»oo poo 

2n(-) m / drr 2 / dt 



oo poo c +iujt 

2 / e 



oo </ — oo 



[ t - (r + ze)] m [ t - (-r + ze)] m ' 
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The t integral in this expression is a contour integral, which is nonzero only for u> > 0, 
where u = w°. For w>0we close the contour upwards in a semicircle, break it into two, 
one around each pole, and use the Cauchy integral formula to get 



1 m {w) = dt 



e +iivt 



[t- (r + ie)] m [t- (-r + ie) 



2ttz Id 171 ' 1 ( e lzuJ \ d m - x ( e iz " \ 

,-rJd^{[z-(r + ie)]m) z _ r+ie > ( } 



(m-1)! I dz™*- 1 \[z - (ie - r)] m J z=r+ie dz 771 ' 1 \[z - (r + ze)]" 



Using the Leibniz rule, d ra /dz ra {f{z) g(z)) = £" =0 {^) f [n ~ s) {z)g {s \z), yields, 

*.(»> = „_ 2 :lS y. ( m : ') (-r <m + » - D! 

p — iruj' 
m+s c 



(ioj) 
x - — - 



m — 1 — s 



(22) 



Now using this in eq. (22) yields, after letting r <-> — r in the second term of (22), 

w- [,j;;; zy E(°:')(-r-i"-^-^ 



s=0 

x^)^- P dr 



(23) 



Finally, since 



r m+s — 2' 
— oo ' 



P/ dr — = 1 J , 
>-oo r (n-1)! 



substitution into (23) yields 

The sum is easily recognised as the second derivative of 2 m ~ 3 [x(l — x)] m_1 , evaluated at 
x = |, and equal to |^r=T ) so we obtain the general expression, 
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So, to summarise, by replacing the delta function by an integral, we are able to do 
each of the k integrals separately, obtaining (18). Collecting the terms in the expansion 
and using (21) and its derivatives to proceed with the x integrals yields the final answer 
to (14). 
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